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= wave  amplitude 

= bottom  area  of  cell 

= cushion  area 

= gap  leakage  area 

= leakage  area 
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= craft  planform  area  (=A  +base  area  of 
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= area  of  orifice  between  loop  and  cell 
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= empirical  coefficients  from  fan  map 
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= depth  of  fully  extended  skirt  below  c.g. 
= cushion  leakage  fraction  of  gap  height 
= heave  displacement 
= moment  of  inertia 
= wavenumber 
= length  of  craft 
= mass  of  craft 
= externally  applied  moment 
= cushion  pressure  (gauge) 

= loop  pressure 
= cell  or  jupe  pressure 
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= volume  rate  of  flow  in,  out  of  cushion 

= length  of  skirt  seal  around  cushion 
peri phery 

= average  skirt  thickness 
= craft  speed 
= unit  step  function 
= cushion  volume 
= craft  weight 
= externally  applied  moment 
- surface  elevation 
= pitch  angle 
= wavelength 
= density 

= angle  of  the  outer  face  of  skirt  with 
respect  to  z-axis 

n 

= encounter  frequency  (=  -y-  UQ+w) 

= wave  frequency 


Subscripts 

b 
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2 


refers  to  bow 

refers  to  the  i**1  seal 

refers  to  jupe  or  cell 

refers  to  plenum  loop 

refers  to  relative  position  with  respect 
to  surface 

refers  to  equilibrium  condition 
(a=0,  h=0,  P=P  = ¥-) 

ftp 

refers  to  forward  subcushion 
refers  to  after  subcushion 
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I.  Introduction 

The  problem  considered  in  this  report  is  the  theoretical  prediction  of  the 
motion  of  an  air  cushion  vehicle  (ACV),  traveling  at  constant  speed  and  course, 
as  it  encounters  head  seas  consisting  of  regular,  sinusoidal  waves.  A mathematical 
model  is  developed  from  first  principles  based  upon  the  dynamic  equations  of  motion 
and  the  air  flow  relations  for  the  air  cushions.  These  equations,  in  their  entirety, 
are  complicated  and  require  numerical  evaluation  by  a digital  computer.  They  have 
thus  formed  the  fundamental  structure  for  computer  simulations  of  ACV  response  to 
waves.  In  this  report,  ho wever,  these  equations  will  be  simplified  and  linearized  so 
that  an  analytical  solution  may  be  obtained.  The  solution  provides  the  pitch  and  heave 
response  of  the  ACV  in  terms  of  the  wave  parameters  and  the  design  particulars  of  the 
craft.  The  purpose  of  this  ideal  model  is  to  determine  the  design  parameters  which 
control  the  various  physical  processes  that  affect  the  seakeeping  ability  of  ACV's. 

The  approach  taken  here  emphasizes  a line  of  research,  in  connection  with  the 
ACV  seakeeping  problem,  which  has  not  been  fully  exploited.  On  the  one  hand,  sophi- 
sticated computer  simulations  are  being  developed  in  which  the  extremely  complex, 
nonlinear  differential  and  algebraic  equations  are  solved  numerically  (1-6).  On  the 
other  hand,  systematic  empirical  studies  are  being  done  utilizing  full  scale  craft 
and  physical  scale  models  in  the  towing  tank  (7-10).  However  relatively  little 
work  has  been  done  using  an  analytical  approach  to  solve  a simplified  set  of  dynamic 
equations.  Though  one  would  expect  that  a closed  form  solution  to  a reduced  set  of 
equations  would  be  less  accurate  than  a numerical  evaluation  of  the  complete,  non- 
linear equations,  the  analytical  approach  has  two  important  benefits.  First  of  all, 
physical  insight  into  the  problem  is  gained  by  close  inspection  and  identification  of 
the  terms  in  the  analysis.  And,  secondly,  recognizable  dynamic  parameters  such  as 
domping  and  restoring  force  terms  can  be  found  explicitly  in  terms  of  the  craft's 
design  particulars.  This,  of  course,  will  be  of  enormous  significance  in  decisions 
made  at  the  design  stage. 
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Efforts  directed  towards  obtaining  analytical  solutions  for  AC V motion  problems 
include  studies  of  ACV  heave  motion  by  Hogben  (11),  Reynolds  (12),  and  Kaplan  and 
Davis  (13).  Lavis,  Bartholomew  and  Jones  (14)  used  a linear  analysis  for  pitch, 
heave  and  roll  motion  in  which  coupled  damped  harmonic  oscillator  equations  were 
employed.  The  critical  stiffness  and  damping  coefficients,  however,  were  specified 
empirically.  Reynolds,  West  and  Brooks  (15)  have  also  developed  a linear  analysis 
for  the  coupled  pitch  and  heave  motion  of  a divided  cushion  ACV  in  head  seas.  In 
this  work,  a rather  limiting  assumption  was  made  that  the  skirt  hem  did  not  come  in 
contact  with  the  wave  surface.  More  recently,  Lebel  and  Swift  (16)  have  obtained 
a solution  to  a nonlinear  model  for  heave  motion,  and  Lundblid  (17)  has  developed  a 
linearized  pitch  and  heave  model  for  single  cushion,  peripheral  all-stabilized  crtft. 

In  t.vis  report,  the  pitch  and  heave  response  of  an  ACV  to  head  seas  comprised 
of  regulai  waves  is  analyzed  using  a theory  which  is  a special  case  of  the  more 
general,  nonlinear  mathematical  models  developed  by  Carrier  and  Swift  (6).  The 
equations  presented  are  a greatly  simplified,  linearized  version  of  the  general  form 
and  are  appropriate  for  instances  in  which  the  dynamic  variables  deviate  only  slightly 
from  their  equilibrium  values. 

Two  craft  are  considered  - a single  cushion,  peripheral  cell-stabilized  ACV 
and  a craft  having  a divided  cushion.  A description  of  the  support  systems  used 
in  each  configuration  is  included  in  reference  (6).  Immediately  following  is  the 
analysis  for  the  single  cushion  c’-aft,  while  the  divided  cushion  ACV  is  considered 
in  the  last  section. 

II.  Single  Cushion,  Peripheral  Cell-Siabilized  ACV 

As  indicated  in  Fig.  1,  the  ACV  is  in  planar  motion  with  its  position  and  orien- 
tation given  by  h and  o,  respectively.  cor  convenience  in  the  analysis,  the  origin 
of  the  x,  2 system,  which  is  fixed  with  respect  to  the  craft,  will  be  at  the  craft 
center  of  gravity  and  over  the  centroid  of  the  planfonu  area.  Whenever  possible, 
advantage  will  be  taken  of  the  smallness  of  the  difference  of  dependent  variables 
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from  their  equilibrium  values  in  order  to  linearize  the  equations.  Products  of 
small  quantities  will  be  dropped  in  favor  of  linear  terms.  The  craft,  traveling 
at  constant  speed  UQ,  encounters  head  seas  consisting  of  regular,  sinusodial  waves 
of  amplitude 


n(x,t)  = a sin(kx  + u t) 


(1) 


The  Froude-Kryloff  hypothesis  will  be  employed;  that  is,  the  waves  are  assumed  non- 
comp li  ant. 

The  equations  of  motion  governing  craft  response  require  that  inertia  forces 
equal  the  forces  on  the  cushion  and  seals.  The  forces  on  the  cushion  or  skirt 
cells  are  taken  as  the  pressure  acting  over  the  externally  supported  area: 


mh  = l [(Pj)k  - PjAb  - [P  - P0]Ap 
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cells 
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skirt 

cells 
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(2) 


(3) 


In  the  above  expressions,  it  has  assumed  that  the  skirt  is  completely  flexible 
and  is  without  inertia.  The  hydrodynamic  pressure  on  the  seal,  when  seal  contact 
occurs,  is  then  equal  to  the  cell  or  jupe  pressure.  With  this  approximation,  the 
effect  of  skirt  contact  is  to  increase  projected  area. 

The  cushion  or  cell  pressures  are  determined  from  a conservation  of  mass  relation. 
For  the  main  cushion,  when  compressibility  effects  are  neglected,  the  rate  of  cushion 
volume  change  must  equal  the  net  flow  into  the  cushion: 
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(4) 


*»  Q.  - Q , 

Min  vout 

The  main  cushion  is  fed  air  directly  by  a fan  supply.  The  linearized  flow-pressure 
relation  represents  the  "fan  map"  near  the  equilibrium  operating  point  and  is  given 
by 


Qin  ‘ Co  " C1P  (5) 

where 


The  flow  out  of  the  main  cushion  is  due  to  leakage  beneath  the  surrounding  skirt. 
An  orifice-coefficient  approach  will  be  used;  in  equation  form  this  is 


«W*V»  (7) 

The  leakage  area,  A£  , includes  a part  present  at  equilibrium  and  additional  area 
due  to  gap  openings  between  skirt  and  surface  which  occur  while  encountering  waves. 

Not  all  the  physical  gap,  however,  is  available  to  main  cushion  leakage.  The  under- 
fed peripheral  jet  flow  from  the  cellular  skirt  must  also  pass  beneath  the  outer  skirt 
along  with  main  cushion  leakage  flow.  Therefore,  only  a fraction  of  the  physical 
gap  area  is  allotted  to  main  cushion  leakage.  The  leakage  area  is  then  given  by 


A£  = Ao  + G v As^pt-h+x^D-f^) 

skirt  f r W 

cells  Vh 


where  G is  the  fraction  gap  area  occupied  by  main  cushion  leakage.  The  fraction  G 
will  be  specified  as  the  ratio  of  leakage  areas,  main  cushion  system  over  skirt 
system  plus  main  cushion  system,  known  to  exist  at  equilibrium.  In  the  design  cen- 
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sidered  the  jupe  pressure  of  the  open-bottomed  cells  at  equilibrium  equals  the  main 
cushion  pressure,  so  that  the  gap  fraction  G may  be  written  as 


G = 


v i <Vo = w; 


(9) 


skirt 

cells 

The  remaining  term  in  Eq.  (4),  the  rate  of  volume  change,  accounts  for  changes  in 
cushion  height  and  the  volume  occupied  by  the  waves  ("wave  pumping").  The  rate  of 
cushion  volume  change  is  then 

xabu  „ 

5 = Ac(-fi) — ■ sin  cosuiet 


(10) 


Combining  Eqs.  (4),  (5),  (7),  (8),  (10),  and  linearizing,  one  obtains 

Aabui  . C A 2P 

M * -r-i  si„  ii  cos»et  = (C,  ♦ )(Mo,  . c„Ag/  01) 

Conservation  of  mass  is  also  applied  to  the  air  flow  through  each  open-bottomed 
cell.  This  is  a statement  that  the  flow  into  the  jupe  from  the  plenum  loop  must  equal 
the  flow  out  through  the  bottom  plus  the  rate  of  volume  change  of  the  cell.  In  equation 
form  this  reads  as 


Vr  — iH—  ■ Ww-ih-  * T 

where  is  the  leakage  area  below  the  outer  skirt  available  to  the  escape  flow, 

Ai  = * (^S5ASi(-h+xi9-n1)Usp(-h+x1-9-ti1) 

and 


02) 


(Ai}o  = V“bjr^ 


(13) 
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Though  the  loop  plenum  Is  fed  separately  from  the  main  cushion,  the  exit  flow 
areas  for  the  cushion  and  the  loop-skirt  system  are  proportional.  The  flow  out  from 
the  two  systems  are  also  nearly  proportional  since  the  equilibrium  pressure  of 
the  open-bottomed  cells  and  the  equilibrium  cushion  pressure  are  the  same.  Because 
of  the  almost  linear  relationship  between  the  flow  through  the  loop-skirt  system 
and  the  main  cushion  system,  the  loop  pressure  will  be  approximated  by 

PL  = V ( 

In  this  expression,  the  constant  of  porportionality,  Cp,  will  equal  the  ratio  of  pres- 
sures at  equilibrium,  Cp  - (PLVPo* 

A L 

The  itn  cell  mass  conservation  relation,  Eq.  12,  is  linearized  by  squaring  both 
sides  and  dropping  terms  involving  the  products  of  the  fluctuations  of  time  dependent 
variables  from  their  equilibriun  values.  After  some  manipulation,  this  may  be  written 
in  the  form 


2(AJ  as. 

<Vfp  ■ - o-ghW-W 


Vo 

(AiV^  „ , , . . , 

m k — T as.  (-n+x^o-n*) 

WW  ’ ’ ’ 

At  the  stern,  the  seals  are  clo'*l  on  the  bottom.  For  these  jupes,  the  cell 
pressure  is  always  equal  to  the  loop  pressure. 

Combining  Eqs,  (11)  and  (15)  with  the  equations  of  motion,  Eqs,  (2)  and  (3), 
one  obtains: 
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(-h+x^o-nj  - PQ  (h-x^e+i^  )U$  (h-x^o+n^)  + 
cot<j>  p 

AK(A.  )n  ^P“ 

+ 2 2 — ^ as-  ("^+x-!®“r'i )} 

cn[A^(Ai  )^3  ’ 1 1 
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(16) 
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It  is  observed  that  these  equations  no*  involve  stiffness  tenas  which  are  piece- 
wise  linear.  That  is,  there  are  now  forcing  terns  of  the  fora 


F<M  “ klhrUsp(+hr^2hrUsp^hri 


{13} 
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where  hr  is  the  relative  skirt  height,  hr  = h-x^e+n^. 

In  order  to  completely  linearize  the  equations,  some  criteria  must  be  set  for  deter- 
mining an  equivalent  stiffness  (spring)  constant  such  that 


F»r>  ■ Keqhr 


(19) 


The  criteria  which  has  given  the  best  results  is  that  the  equivalent  stiffness  be 
such  that  the  period  of  free  oscillation  be  the  same  for  the  Eq.  (19)  form  as  for 
the  Eq.  (18)  form.  This  criteria  implies  that 


K 


4k-j  k2 


‘eq 


k-j+2  /k-|  kg+kg 


(20) 


After  this  approximation  is  made,  the  summations  may  be  evaluated  by  calculating 
the  corresponding  integrals  along  the  craft  periphery.  The  equations  can  then  be  put 
in  the  linearized  form 


(m+A33^+B33^+C33h+A356+B359+C3:>e  ' Zext 


A53^+B53^+C53h+^y+A55^®+B55®+C559  = Mext. 
in  which  the  hydrodynamic  coefficients  and  wave  forcing  terms  are  given  explicitly 

in  Appendix  A. 

Since  the  wave  forcing  is  a harmonic  function  of  time  at  the  encounter  fre- 
quency, Zext  and  Mext  can  be  expressed  in  the  form 


(21) 


lVt 


Zext  = Re^Zoe  1 


iw  t 

Hext  “ W 6 3 


(22) 


where  ZQ  and  MQ  are  the  complex  amplitudes  of  the  wave  force  and  moment  on  the  craft. 


The  craft  response,  h(t)  and  e(t),  are  expressed  similarly: 

iiu.t 


h(t)  - Re  [K0e  • ] 


1(0  t 

e(t)  = Re  [eQe  ] 


(23) 


n 


Substituting  Eqs.  (22)  and  (23)  into  Equation  (21)  enables  the  complex  amplitudes 
of  the  response,  hQ  and  8 , to  be  determined  in  the  frequency  domain.  The  results 


are: 


ho  = {Zo  [C55"“e(Iy+A55)+i(1)e  B55]"Mo[C35'(i)eA35+ia)B35]}* 

‘^C33‘we^m+A33^  + iue  B33^C55"“e^y+A55^  + iu)e  B55^ 

"^■C53““e  A53  + 1ai  B53^C53"ue  A35  + iwe  B35^ 

eo  = {^o^C33"“e^m+A33^  + ia)e  B33^“^o^C53"“e  A53+1a>e  B53^‘ 

•{[C33-aje  (m+A33)  + i^e  B33][C55  - (Iy+A55)  + i^e  B55] 


'^C53  " we  A53  + iwe  B53^C35  " “e  A35  + iwe  B35^ 


(24) 


Eq.  (23)  and  (24)  rer”esent  the  solution  for  the  craft  response  variables  e(t)  and  h(t). 
These  expressions  ar«  written  in  terms  of  dynamic  parameters  which  are  given  explicitly 
in  Appendix  A as  functions  of  design  particulars. 


III.  Divided  Cushion  ACV 

The  pitch  and  heave  motion  of  a divided  cushion  ACV  will  be  analyzed  using 
the  same  coordinate  system  as  that  used  in  the  previous  section  (i.e..  Fig.  1).  The 
vehicle  considered  here  will  have  its  rectangular  cushion  divided  by  a transverse 
stability  keel  at  a point  midway  between  the  bow  and  stern  coinciding  with  the  craft 
center  of  gravity.  (A  longitudinal  stability  keel  has  no  significance  in  a head 
seas,  pitch  and  heave  analysis).  The  forward  cushion,  designated  with  subscript  1, 
and  the  after  cushion,  designated  by  subscript  2,  are  both  enclosed  externally  by  a 
flexible  skirt  consisting  of  finger  seals.  Thus,  there  is  no  cell  base  area  and 


The  craft,  traveling  at  constant  speed  UQ,  encounters  the  regular  wave  head  sea 
given  by  Eq.  (1).  Assuming  that  the  skirt  system  is  completely  flexible  and  without 
inertia,  the  equations  of  motion  for  pitch  and  heave  can  be  written  in  the  form: 
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«*  = - (p]-po>  Ac/2  - (p2'p0)  Ac/2  • J (h-Xj e+tij )U$p (h-Xi e+n, ) 4s, 


external 

periphery 


y ■ C(prpo>  ' (p2-po>^+  \ ^xl(h-xie+pf)usp<h-x1e+''i> 


(25) 


external 

periphery 


The  pressure  fluctuations  from  their  equilibrium  values  are  determined  by  con- 
sidering the  flow  relations  for  each  cushion.  Since  the  leakage  areas  are  large,  the 
density  is  assumed  constant,  and  the  rate  of  cushion  volume  change  equals  the  net 
flow  into  the  cushion: 


V2  = (Qinh,2  ’ ((Wl,2 


(26) 


In  the  craft  considered,  each  cushion  is  supplied  air  by  a separate  fan  system.  The 
linearized  representation  of  the  fan  map  is  then  given  by 


^1nh,2  ’ Co  " C1  pl,2 


(27) 


The  flow  out  is  expressed  as  an  orifice  coefficent  modification  to  Bernoulli's 
equation: 


2T 


((W\2  = Wl,2/ — ^ 


Wl,2/T  ( } 

° 2^ 


where 


(A. 


VC1  po 


ih.2  “ 1WC  * 1 asi  <*hr,1  usp  f'Vi 

r / C 


c / 


(28) 


n p external 

periphe»*y  1,2 


In  the  above  equation,  cross-flow  leakage  between  the  cushions  is  assumed  negllbible 
in  comparison  to  leakage  to  the  atmosphere. 
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The  rate  of  cushion  volume  change  is  due  to  changes  in  height  of  the  hardstructure 
base  minus  the  rate  of  volume  change  occupied  by  the  wave.  In  equation  form  this  is 

xabu)  n « 

^ 2 = (“f*  ± £/4  ®)(Ac/2) “ sin  fx  cos  (ff  1 “e1^  (29) 

Combining  Eqs.  (26)-(29),  one  obtains  the  pressure  fluctuations  from  equilibrium 


as 


P,  , - Pn  = 

1,2  0 


1 


r j.  C -C,Pn 

Ci+  0 1 0 


{(fm/4e)  Ac/ 2 + 


Xabo), 


sinf  cos(f  1Wet) 


2F 


rr 

- c ✓ — - 

n p 


l Asi (“hr ) iUSp(~hr) i ) 


external 
periphery  1,2 


(30) 


Substituting  this  result  into  the  equations  of  motion,  Eq.  (25),  yields, 


mh 
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0 1 0 
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external 
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(31) 


It  *s  seen  that  the  restoring  force  terms  are  piecewise  linear,  the  general  fora 


of  which  is  given  by  Eq.  (18).  For  the  divided  cushion  craft,  a generating  function 
method  was  used  to  determine  the  equivalent  stiffness  constant  defined  by  Eq.  (19). 
This  criteria  requires  that  an  average  of  the  stiffness  constants  be  taken: 


Kj  + Kg 


(32) 


By  using  this  approximation  ana  evaluating  the  summations  by  determining  the  corres- 
ponding integrals,  Eq.  (31)  is  put  in  the  form  of  Eq.  (21)  in  which  the  coefficients 
are  given  explicitly  in  Appendix  B.  The  solutions  for  the  pitch  and  heave  variables, 
e(t)  and  h(t)  are  then  given  by  Eqs.  (23)  and  (24). 
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APPENDIX  A 

The  coefficients  and  wave  forcing  terms  in  Eq.  (21)  for  the  single  cushion, 
peripheral  cell-stabilized  craft  are  as  follows: 
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APPENDIX  B 


The  coefficients  and  wave  forcing  terms  in  Eq.  (21)  for  the  divided  cushion 
ACV  are  as  follows: 
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